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Abstract. The paper presents a methodology to enhance the stiffness analysis of parallel manipu- 
lators with parallelogram-based linkage. It directly takes into account the influence of the external 
loading and allows computing both the non-linear "load-deflection" relation and relevant rank- 
deficient stiffness matrix. An equivalent bar-type pseudo-rigid model is also proposed to describe 
the parallelogram stiffness by means of five mutually coupled virtual springs. The contributions 
of this paper are highlighted with a parallelogram-type linkage used in a manipulator from the 
Orthoglide family. 
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ing, linear approximation. 



1 Introduction 

In the last decades, parallel manipulators have received growing attention in indus- 
trial robotics due to their inherent advantages of providing better accuracy, lower 
mass/inertia properties, and higher structural rigidity compared to their serial coun- 
terparts (Ceccarelli et al., 2002, Company et al., 2002, Merlet, 2000). These features 
are induced by the specific kinematic structure, which eliminates the cantilever-type 
loading and decreases deflections due to external wrehches. Accordingly, they are 
used in innovative robotic systems, but practical utilization of the potential benefits 
requires development of efficient stiffness modeling techniques, which satisfy the 
computational speed and accuracy requirements of relevant design procedures. 

Amongst numerous parallel architectures studied in robotics literature, of special 
interest are the parallelogram-based manipulators that employ special type of link- 
age constraining undesirable motions of the end-platform. However, relevant stiff- 
ness analysis is quite complicated due to the presence of internal closed kinematic 
chains that are usually replaced by equivalent limbs (Majou et al., 2007, Pashkevich 
et al., 2009a) whose parameters are evaluated rather intuitively. Thus, the problem 
of adequate stiffness modelling of parallelogram-type linkages, which is in the focus 
of this paper, is still a challenge and requires some developments. 
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Another important research issue is related to the influence of the external load- 
ing that may change the stiffness properties of the manipulator In this case, in ad- 
dition to the conventional "elastic stiffness" in the joints, it is necessary to take into 
account the "geometrical stiffness" due to the change in the manipulator configura- 
tion under the load (Kovecses and Anjeles, 2007). Moreover buckling phenomena 
may appear (Timoshenko and Goodier, 1970) and produce structural failures, which 
must be detected by relevant stiffness models. 

This paper is based on our previous work on the stiffness analysis of over- 
constrained parallel architectures (Pashkevich et al., 2009a, 2009b) and presents 
new results by considering the loading influence on the manipulator configura- 
tion and, consequently, on its Jacobian and Hessian. It implements the virtual joint 
method (VJM) introduce by Salisbary (Salisbury, 1980) and Gosselin (Gosselin, 
1990) that describes the manipulator element compliance with a set of localized 
six-dimensional springs separated by rigid links and perfect joints. The main con- 
tribution of the paper is the introduction of a non-linear stiffness model of the 
parallelogram-type linkage and its linear approximation by 6x6 matrix of the rank 
5, for which it is derived an analytical expression. 



2 Problem statement 



Let us consider a typical parallel manipulator that is composed of several kinematic 
chains connecting a fixed base and a moving platform (Figured^). It is assumed, that 
at least one chain includes a parallelogram-based linkage that may introduce some 
redundant constraints improving the mechanism stiffness. Following the VJM con- 
cept, the manipulator chains are usually presented as a serial sequence of pseudo- 
rigid links separated by rotational or translational joints of one of the following 
types: (i) perfect passive joints; (ii) perfect actuated joints, and (iii) virtual flexible 
joints that describe compliance of the links, joints and/or actuators. 




To derive a parallel architecture without internal kinematic loops, the parallel- 
ogram based linkage must be replaced by an equivalent bar-type element that has 
similar elastic properties. Thus, let us consider the VJM model of the kinematic par- 
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allelogram assuming its compliance is mainly due to the compliance of the longest 
links, which are oriented in the direction of the linkage (Figure[T]3,c). In order to be 
precise, let us also assume that the stiffness of the original bar-elements is described 
by a 6 X 6 matrix whose elements are identified using FEA-based modeling. 

In the framework of the VJM approach, the geometry of the parallelogram-based 
linkage can be descried by the following homogeneous transformations 

T, = T,(77 d/2) R.fei) T,(L) T,(0,i) Tv(0,2) T,(0,-3) ^M^) 

R,(0,-5) Rz(0/6) R,fe2) T,(-T7 d/l) K,{-qa) 

where r\ = (—1)', T, is a 4 x 4 homogenous transformation matrix, T^, R,. are 
the matrices of elementary translation and rotation, L, d are the parallelogram ge- 
ometrical parameters, qtj are the passive joint coordinates, 0,-, are the virtual joint 
coordinates, i — 1,2 defines the number of the kinematic chain, j identifies the co- 
ordinate number within the chain. It should be noted that using this notation, the 
closed loop equation can be expressed as Ti = T2 . 

For further computational convenience, the homogenous matrix equations (1) 
may be also rewritten in the vector form as 

t = g/(q„ft), (2) 

where vector t = (p, <p) is the output frame pose, that includes its position p — 
(x^y^zY and orientation <p = {(p^, (Py, (Pz), vector q/ = contains the passive 

joint coordinates, and vector 0,- ~ {On , Q/g) collects all virtual joint coordinates of 
the corresponding chain. Using the above assumptions and definitions, let us derive 
the stiffness model of the parallelogram-based linkage. 



3 Static equilibrium 

Let us derive first the general static equilibrium equations of the parallelogram as- 
suming that it is virtually divided into two symmetrical serial kinematic chains with 
the geometrical and elastostatic models t = g,(q,, 0;), T, = K,- 0,-, / = 1,2 where the 
variables 0, correspond to the deformations in the virtual springs, K,- is the spring 
stiffness matrix, and t,- incorporates corresponding forces and torques. Taking into 
account that the considered mechanism is under-constrained, it is prudent to assume 
that the end-points of both chains are located at a given position t and to compute 
partial forces F, corresponding to the static equilibrium. 

To derive equations, let us express the potential energy of the mechanism as 



^ 1=1 



(3) 
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In the equilibrium configuration, this energy must be minimised subject to the geo- 
metrical constraints t = g,(q,; ft), i — 1,2. Hence, the Lagrange function is 

1 2 2 

, 02,qi ,q2) = 2 L ^ei ft + L (t - gi{0„qd) (4) 

/= 1 /= 1 

where the multipliers A, may be interpreted as the external forces F,- applied at the 
end-points of the chains. Further, after computing the partial derivatives of L{...) 
with respect to ft, q,-,A,- and setting the derivatives to zero, equations of the static 
equilibrium can be presented as 

J^,.A,- = K; ft; 4 A, = 0; t = g,(q,-, ft); /=1,2 (5) 

where Jg,- — (9g,(...)/(9ft and J^,- = 5g,(...)/(9q, are the kinematic Jacobians. Here, 
the configuration variables ft-, q, and the multipUers A,- are treated as unknowns. 

Since the derived system is highly nonlinear, the desired solution can be obtained 
only numerically. In this paper, it is proposed to use the following iterative procedure 



q'i 







t-g/ + J</iq; + Je; ft 




0/ = Kfl/jS,A',- (6) 



where the symbol " ' " corresponds to the next iteration. Using this iterative pro- 
cedure, for any given location of the end-point t, one can compute both the partial 
forces F, and the total external force, which allows to obtain the desired "force- 
deflection" relation F — f (t) for any initial unloaded configuration. 



4 Stiffness matrix 

To compute the stiffness matrix of the considered parallelogram-based mechanism, 
the obtained "force-deflection" relations must be linearized in the neighborhood of 
the static equilibrium. Let assume that the external forces F,, / = 1,2 and the end- 
point location t of both kinematic chains are incremented by some small values 
5F,, 8i and consider simultaneously two equilibriums corresponding to the state 
variables (F, , ft ,q,-,t) and (F, + 5F,, ft + 5ft ,q, + 5q, ,t + 5t). Under these assump- 
tions, the original system (|5]) should be supplemented by the equations 

(Je; + dieif (A,- + 5A,) = Kg,- ( ft + 5ft) ; 

{igi + 5i,i) (Ai + 5A,-) = 0; (7) 

t + 5t = g,(q,-, ft) + Je; 5ft + 5q,-; 

where 5Je,- and 51^, are the differentials of the Jacobians due to changes in (ft,q,). 
After relevant transformation, neglecting high-order small terms and expanding the 
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differentials via the Hessians of the function = gi(qj, 0,)^ A, ; H^'^ — d^^i/dqj, 
Hg = d^Wi/dqfOi, = d^Wi/ddf, the system of equations may be rewritten as 



limination of 
one can obtain a matrix equation 



50, = 0; 



(8) 



After analytical elimination of the variable 5 0, and defining k^^ = (Kg — Hg^ ) ' , 



"5A,' 




5q, 





Je; • k' 



(0 



iqi + Je; ■ kg"* 



H 



(0 

9q 



which yields the linear relation 5A = K^. 



-, -1 


"5t' 








(9) 



5t defining the Cartesian stiffness matrix 



Kc, for each kinematic chain. Taking into account the architecture of the consid- 
ered mechanism, the total stiffness matrix may be found as = K^i + ^c2- These 
expressions allow numerical computation of the Cartesian stiffness matrix for a gen- 
eral case, both for loaded and unloaded equilibrium configuration. 

However, in the case of unloaded equilibrium, the above presented equations may 
be essentially simphfied: 
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(10) 



where Cq = cosq, Sq = sinq, S2q — sin 2^, q = qu (Figure lb,c), Kij are the elements 
of the 6x6 stiffness matrix of the parallelogram bars K^, which are assumed here 
to be the only elements of the linkage that posses non-neghgible stiffness. In the 
following example, this expression is evaluated from point of view of accuracy and 
applicabihty to stiffness analysis of the parallelogram-based manipulators. 



5 Illustrative example 

To demonstrate validity of the proposed model and to evaluate its applicability 
range, let us apply it to the stiffness analysis of the Orthoghde manipulator (Fig- 
ure la). It includes three parallelogram-type linkage where the main flexibility 
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source is concentrated in the bar elements of length 310 mm. Using the FEA-based 
software tools and dedicated identification technique (Pashkevich et al., 2009c), the 
stiffness matrix of the bar element was evaluated as 



K,, 



"2.20- 10-* 




















1.81 ■ IQi 











-2.84 -10^ 








7.86-10' 





1.25-10'* 
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-2.84-103 











5.85-10' 



(11) 



where for linear/angular displacements and for the force/torque are used the fol- 
lowing units: [mm], [rad], [N], [Nmm]. This bar-element was also evaluated with 
respect to the structural stability under compression, and the FEA-modelling pro- 
duced three possible buckling configurations presented in Figure 2. It is obvious 
that these configurations are potentially dangerous for the compressed parallelo- 
gram. However, the parallelogram may also produce some other types of buckling 
because of presents of the passive joints. 




Fig. 2 Critical forces and buckling configurations of a bar element 
employed in the parallelograms of Orthoglide. 



For the parallelogram-base linkage incorporating the above bar elements, expres- 
sion ( fTTT i yields the following stiffness matrix 
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5.92 -10' 



(12) 



corresponding to the straight configuration of the linkage (i.e. to ^ = 0). Being in 
good agreement with physical sense, this matrix is rank deficient and incorporates 
exactly one zero row/column corresponding to the z-translation, where the parallelo- 
gram is completely non-resistant due to specific arrangement of passive joints. Also, 
as follows from comparison with doubled stiffness matrix ( fTTI ) that may be used as a 
reference, the parallelogram demonstrates essentially higher rotational stiffness that 
mainly depends on the translational stiffness parameters of the bar (moreover, the 
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rotational stiffness parameter ^^55 of the bar is completely eliminated by the pas- 
sive joints). The most significant here is the parallelogram width d that explicitly 
presented in the rotational sub-block of Kp. 

To investigate applicability range of the linear model based on the stiffness matrix 
([Tol l, it was computed a non-linear "force-deflection" relation corresponding to the 
parallelogram compression in the x-direction. This computation was performed us- 
ing an iterative algorithm presented in Section 3. As follows from obtained results, 
the matrix (fTol i ensures rather accurate description of the parallelogram stiffness 
in small-deflection area. However, for large deflections, corresponding VJM-model 
detects a geometrical buckling that limits applicability of the matrix ( fTOl i. 

Similar analysis was performed using the FEA-technique, which yielded almost 
the same "force-deflection" plot for the small deflections but detected several types 
of the buckling, with the critical forces may be both lower and higher then in the 
VJM-modelling. Corresponding numerical values and parallelogram configurations 
are presented in Figure 3. 

The developed VJM-model of parallelogram was verified in the frame of the stiff- 
ness modelling of the entire manipulator Relevant results are presented in Table 1 . 
They confirm advantages of the parallelogram-based architectures with respect to 
the translational stiffness and perfectly match to the values obtained from FEA- 
method. However, releasing some assumptions concerning the stiffness properties 
of the remaining parallelogram elements (other than bars) modifies the values of 




Fig. 3 Force-deflection relations and buckling configurations for parallelogram compression 
(modelling methods: VJM, FEA, and LIN - linear model with proposed stiffness matrix). 



Table 1 Stiffness parameters of the Orthoglide manipulator for different assumptions concerning 
parallelogram linkage. 



Model of linkage 




VJM-model 




FEA-model 






Kmt 




Kmt 




[mm/N] 


[rad/N] 


[mm/N] 


[rad/N] 


2x-bar linkage 


3.35- 10^ 


0.13-10*' 






Parallelogram linkage 










(solid axis) 


3.23- 10^ 


4.33 - 10*' 


3.33- 10^ 


4.10 - 10*5 


Parallelogram linkage 










(flexible axis) 


3.08-103 


4.06-10*^ 


3.31 - 10^ 


4.07 - 10*^ 
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the translational and rotational stiffness. The latter gives a new prospective research 
direction that targeted at more general stiffness model of the parallelogram. 



6 Conclusions 

The paper presents new results in the area of enhanced stiffness modeling of par- 
allel manipulators with parallelogram-based linkage. In contrast to other works, it 
expUcitly takes into account influence of the loading and allows both computing 
the non-linear load-deflection relation and detecting the buckling phenomena that 
may lead to manipulator structural failure. There is also derived an analytical ex- 
pression for relevant rank-deficient stiffness matrix and proposed equivalent bar- 
type pseudo-rigid model that describes the parallelogram stiffness by five localized 
mutually coupled virtual springs. These results are vahdated for a case study that 
deals with stiffness modeling of a parallel manipulator of the Orthoglide family, for 
which the parallelogram-type linkage was evaluated using both proposed technique 
and straightforward FEA-modeling. 
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